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Abstract 

The first linear global electromagnetic gyrokinetic particle simulation on the excita- 
tion of toroidicity induced Alfven eigenmode (TAE) by energetic particles is reported. 
With an increase in the energetic particle pressure, the TAE frequency moves down 
into the lower continuum. 

PACS numbers: 52.55.Fa, 52.35.Bj, 52.65.Tt 

Toroidicity induced Alfven eigenmode (TAE)^^^ can play important roles in burning 
plasmas. The TAE modes can be excited when energetic particles, for example fusion born 
alpha particles, resonate with the phase velocity of the shear Alfven wave which resides 
within the frequency gap of the Alfven continuum. 

Shear Alfven wave oscillations, continuum damping, and the appearance of the frequency 
gap in toroidal geometries by gyrokinetic particle simulation have been recently reported.^ 
The simulation of Ref.4 is demonstrated in the long wavelength magnetohydrodynamic 
(MHD) like limit in the absence of kinetic ions. In this letter, taking exactly the same 
parameters''^ but adding the energetic ion particles, the first linear particle simulation on 
the excitation of the TAE modes is reported. The simulation is done without employing 
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MHD equation. The simulation is not the conventional gyrokinetic-MHD hybrid ones,^~^ 
where the kinetic ions enter the system through the stress tensor. The setting of the sim- 
ulation is kept as faithful as possible to Ref.3 and 4 to see an explicit connection with our 
previous studies.^ 

A simplified linearized set of equations is employed^ for the numerical simulation which 
is reduced from the electron-fluid ion-kinetic hybrid gyrokinetic model. ^^^^ The equations of 
Ref.4 are normalized by the ion Larmor radius (at the electron temperature) for the length, 
the ion cyclotron frequency for time, and the electron temperature for the electrostatic 
potential, and the magnetic field strength at the magnetic axis, Bq. The set of the equations 
are the electron continuity equation 

^ = -V„fc,|. (1) 

{Srie is the fluid electron density and Su^e is parallel electron velocity), the inverse of Fara- 
day's law 



= V|| ($e// - ^) (2) 



dt 

{A^^ is the vector potential, $ is the electrostatic potential, and $e// is the effective potential 
representing the total parallel electric field) , the gyrokinetic Poisson equation^^ 

$ - $ = Suc, - 6ne (3) 

{Sria is the gyro-averaged energetic particle density, $ is the second gyrophase averaged 
electrostatic potentiaP^) the lowest order adiabatic relation 

^eff = Srie, (4) 
and the inverse of Ampere's law 

6uile = /3e"^Vi^|| + 6u\\a- (5) 

The parallel velocity of the energetic particles is given by Su\\a. Here f3e = (cs/va)'^ where Cg 
is the sound velocity and va is the Alfven velocity. All the variables in Eqs.(l)-(5) are the 



normalized ones. The gradient operators V_l and Vy are those in the direction perpendicular 
and parallel to the equilibrium magnetic field. 

By coupling Eqs.(l)-(5), the shear Alfven wave dispersion relation in the toroidal geom- 
etry, Eq.(2) of Ref.3 can be obtained. Figure 1 shows the shear Alfven wave frequency as 
a function of the radial coordinate r (a is the minor radius), which is equivalent to Fig.l 
of Ref.3. Due to the 1/ R variation of the toroidal magnetic field {R is the major radius), 
the cylindrical Alfven continuum (dashed lines) breaks up and the frequency gap (or the 
frequency forbidden band) appears within the range of 0.299 < uj/uja < 0.389. Here, ua is 
the Alfven frequency at the magnetic axis. 

Equations (l)-(5) are employed (with the Sua and Su\\a terms turned ojf) in the simu- 
lation of Fig. 5 of Ref. 4. On top of Eqs.(l)-(5) we add kinetic ions. The guiding center 
equation and the 5f gyrokinetic equation (the weight equation) are solved for the kinetic 
energetic particle ions (we neglect thermal kinetic ions, however). Taking the perturbed 
distribution function 5 fa, the energetic particle density (velocity) in Eq.3 [Eq.(5)] is given 
by = TaJ 6fad^v (5M||a = Ta J V\\6fad^v), where / d^v is an average over the velocity 
space. Note that in the simulation, we control the perturbed density of the energetic parti- 
cles by multiplying a factor r„ (r^ < 1) which is proportional to the equilibrium energetic 
particle density. The energetic ion particles are provided with the Maxwellian distribution 
function /qq oc exp (— Vy/2v^) in the velocity space (thus dfoa/dv\\ is always negative). The 
resonating energetic particles are inserted within the frequency gap by utilizing a portion of 
the Maxwellian distribution function whose thermal velocity is of the order of Alfven speed. 

The particles that resonate with the shear Alfven wave with the phase velocity cu/k^^ 
can destabilize the TAE mode, when the mode frequency cu is within the frequency gap (we 
choose uj/uja = 0-344 in the middle of the gap in Fig.l), and when the parallel wave vector 
/cy — {m — nq)/qR satisfies k\\ — —k\\m — ^||m+i at g = (2m -|- l)/2n. Here, m {n) stands 
for the poloidal (toroidal) mode number. We take m — 1, m -\- 1 — 2, and n — 1 which is 
equivalent to m = —2, m + 1 = —1, and n = — 1 of Ref.3. The geometrical parameters used 
for the simulation are the same as in Refs.3 and 4 (for example, the inverse aspect ratio of 
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0.375 and a parabolic safety factor q). The major radius is given by it! = 46.6cm as well 
(after convincing the TAE excitation in the originally published setting,^ we move on to a 
parameter survey in a larger size plasma). From the cu and the k\\ values chosen above, we 
provide the Maxwellian distribution with Va — = 10.32c5. The mass and the charge 
of the energetic particles are that of the Hydrogen ion. In the specific simulation below, 
we set = 0.01, and the constant density gradient «;„ = —R{l/na){dna/dr) = 8.0. Here, 
Ua represents the equilibrium density of the energetic particles. The temperature gradient 
parameters^ are set to be zero. In Eq.(3) and Eq.(5), = 0.15 is taken for Figs. 2 and 3. 

The simulation is conducted by an electromagnetic extension*^ of the GTC code^^^^^ with 
a non-iterative field solver. With the additional energetic particle drive, the TAE mode 
is excited. A linear eigenmode (contour plot) of the TAE instabihty is shown in Fig.2. Note 
that the contour plots are not up-down symmetric due to the finite poloidal shear fiow.^^ 

The frequency spectrum of the TAE instabihty is shown in Fig. 3. The TAE frequency 
{u/loa — —0.36) is found within the gap (and not on the gap boundaries as in Ref.4)^^ which 
is a clear evidence of the TAE excitation. The linear growth rate of the TAE instability is 
given by ^/coa — 0.0215 (and thus \^/u;\ — 6.0%) for both the m — 1 and m — 2 mode. 

Figure 4 shows the hnear TAE growth rates (divided by the real frequency) versus the 
multiplication factor r^. Compared to the calculations in Figs. 2 and 3, a twice larger plasma 
size is taken for Fig.4 (the Larmor radius of the energetic particles in Figs. 2 and 3 are 15% 
of the minor radius while 7.5% in Fig.4). We see a monotonic increase in the growth rate as 
the energetic particle population increase (and thus the effective beta value of the energetic 
particles, (3a increases; at Pe — 0.01, a simple estimate will give Pa — '^T^'^aTa/ Bq ~ Tq,. Here 
Ta is the energetic particle temperature). On the other hand, the real frequency of the mode 
decreases (approximately 15% of a reduction in the real frequency) as (or Pa) increases 
and crosses the lower gap boundary^^ (but not the upper gap boundary) which is suggested 
by the analysis in Ref.24 (Resonant TAE, which emphasizes the resonance between the mode 
frequency and the magnetic drift frequency of the energetic particles). As a reference, the 
energetic particle mode (EPM)^^ refers to a heating of the continuum based on the notion 
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that the energetic particle drive exceeds the continuum damping and predicts the appearance 
of the mode frequency both in the upper and the lower continuum. The square plots in Fig. 4 
represents the analysis of the TAE growth rate in a large aspect ratio tokamak from Ref.3 
(the simulation results and the analysis compare favorably at higher (5q). 

We note that instability growth was already minimal at — 0.025 (with the specific 
simulation parameters we employed in Fig. 4), and we did not survey below < 0.025 in this 
work. The TAE mode in its nature should not have an instabihty threshold. The latter onset 
feature needs to be investigated in detail to see the limitation of the initial value approach 
(if any). We also would like to remind that a simplified model Eqs.(l)-(5) is employed in 
this letter (so as to primary focus on the excitation of TAE by the additional energetic 
particles).^ The radial extension of the simulation domain is limited to 0.1 < r/a < 0.9 (see 
Fig. 2). An inclusion of the magnetic axis can be crucial to describe the long wavelength 
global modes precisely. 

In summary, the first linear excitation of the low-n TAE modes by the energetic particles 
in a global gyrokinetic particle simulation is reported. The work did not employ MHD 
model (through closure relations). With a completion of the current global gyrokinetic 
simulation method, one can investigate the onset and the saturation mechanism of the TAE 
modes simultaneously without any restrictions on the wavelength of the modes. Apparently, 
the advantage of initial value approach is its application for nonlinear simulation. We plan to 
report the analysis of energetic particles driven high-n Alfvenic modes separately. Whether 
which mode numbers are most unstable is a great interest to large tokamak burning plasma 
experiments. 
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FIG. 1. Shear Alfven frequency as a function of the radial location. The dashed curves are the 
continuum frequencies for the cylindrical limit, for m = 1 and m = 2 modes. The solid lines are 
for the continuum frequency with the toroidal geometry effect. The lower (upper) boundary of the 
upper (lower) curve is at u/loa = 0.389 {lo/ua = 0.299). Correspondingly the frequency gap (the 
forbidden frequency range) appears within the range of 0.299 < lo/uja < 0.389. The figure is a 
recapitulation of Fig.l of Ref.3. 
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FIG. 2. (Color online) Linear eigenmodes (contour plots on a poloidal plane) of TAE instability, 
(a) The electrostatic potential $ at a toroidal angle C = 0- R^d represents positive $ values, while 
blue represents negative $ values, (b) The vector potential A\\ at a toroidal angle C = 0- 
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FIG. 3. (Color online) (a) The Fourier components of the electrostatic potential ^ as a function 
of time, (b) The frequency spectrum obtained from the time series of Fig. 3(a) (black for m = 1, 
red for m = 2). 
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FIG. 4. (Color online) Dependence of the real frequency w, and the linear growth rate 7 on 
(and thus Pa)- For each data, eigenmode structures similar to Fig. 2 are obtained. The green 
squares represent the analytical TAE growth rate - real frequency ratio in a large aspect ratio 
tokamak from Ref.3. 
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